Vanishing dynamical quark mass at zero virtuality? 
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We show that the dynamical quark mass in effective nonlocal models can van- 
ish at zero virtuality of the quark as M(p^) oc p^. Our arguments follow from the 
constrained-instanton model of the QCD vacuum and from QCD sum rules calcu- 
li , lations with nonlocal condensates. The discussed models also lead to analyticity of 
I M(p^) in the vicinity of zero. 
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The concept of dynamical quark mass generated by the spontaneous breaking 
of chiral symmetry has proven very useful in the description of low-energy phe- 
nomena. This mass enters the quark propagator, Sp-^i^p) = Z{jP')[^ + M(p^)], 
and, in general, depends on the (Euclidean) momentum p. Specific forms of 
M{p^) are obtained in instanton models |T|] or via the Schwinger- Dyson resum- 
mation techniques [Q]. The quark mass in its own is not observable and does 
not bear physical significance. Moreover, it depends on the choice of the gluon 
field gauge, hence is not unique. Also, the region where, e.g., the instanton 
models dominate is at some intermediate values of p^, not too high, where the 
model predicts a too fast dropping of the mass in the perturbative region, and 
not too low, since certainly the infrared physics (confinement) is not being 
incorporated properly. Nevertheless, the mentioned popular approaches have 
founded a typical image of M{p^) in the Euclidean space: it is non-zero at 



* Research supported in part by the Polish State Committee for Scientific Research 
grant 2P03B-080-12, the Russian Foundation of Basic Research grant 01-02-16431, 
INTAS-00-00366, and by the Bogoliubov-Infeld programm. 

E-mail: dorokhov@thsunl.jinr.ru, broniows@solaris.ifj.edu.pl 



= 0, drops monotonically with increasing p^, and asymptotically reaches 
the current quark mass value, m^.Q 



In this paper we argue that this conventional picture need not be the case, 
and that we may well have 

M(p2 = 0) = 0. (1) 

As is increased, M(p^) grows, arrives at maximum, and then decreases to 
reach rric asymptotically. In fact, we will present arguments that the situation 
(J^) can be more natural than the conventional case of M(p^ = 0) > 0. One 
class of arguments comes from QCD sum rules with nonlocal condensates, an- 
other one from the so-called "constrained" instantons, where the large-distance 
asymptotics is improved by an additional cut-off function in the instanton size. 

We will also show that although Eq. (P looks, at first glance, unusual, it 
does not lead to pathologies: the chiral symmetry is broken, and all results 
enforced by this fact hold: the quark condensate has a non-zero value, the 
Gell-Mann-Oaks-Renner relation is satisfied, as are other relations coming 
from symmetries, anomalies are preserved, etc. 

We begin by presenting a calculation in framework of the constrained instanton 
model, where Eq. (|l]) holds. We use for simplicity the "constrained" quark zero 
mode in the singular gauge and introduce the following ansatz for the profile 
function 0]: 

where x is a color Dirac spinor given by X^X^ = (7/i7i^/16) (1 ± 75) /2r^r^ 
and = (=F^, with the upper (lower) signs corresponding to solutions 
in the instanton (anti-instanton) field. The standard instanton solution is ob- 
tained with the constant for p^(x^) = p^. For the constrained instantons one 
uses exponentially-decreasing functions fix), normalized as p^(0) = p^. We 
shell use the form [ffl 
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where K^, (z) is the modified Bessel function. The specific feature of the con- 
strained instanton is that at small distances it is close to the standard in- 
stanton profile of size p, and at large distances it has exponentially-decreasing 
asymptotics governed by a large-scale parameter R, such as p < R. These 
shapes are motivated by considering modifications of the instanton in the 



Throughout this paper we are working in the strict chiral limit of rric = 0. 
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background field of large-scale vacuum fluctuations [Q. The constrained in- 
stanton profile, as opposed to the unconstrained one, provides the correct 
large-distance asymptotics of the quark and gluon field correlators [Q. 

Let us define the 4-dimensional Fourier transform of the quark zero mode 
profile <fsing{p)- The explicit form of the Fourier transform of the standard 
(unconstrained) zero mode in the singular gauge is well known, 

^LM) = [hiz)K^iz) - Hz)K,{z)U^„ , (4) 

and has the asymptotics 

1 ^ as p2^0, 
^\p-) = { ; \ (5) 

^ as p' ^ oo. 

The constrained zero mode has the same large-p^ asymptotics, but it goes to 
a constant at _ > Q 

ip^\p-) = { J"' I (6) 

^ as p^ oo, 

where Ici = Sq" du uipciiu). The quark mass generated by the one-instanton 
background is equal to p 

M(/)=CpV'(p'), (7) 
where the constant C > is determined from the gap equation 

d^k M'^ik'^) n 



(27r)4p + M2(A;2) 4N^' ' ' 

with n = 0.0016 GeV^ ^ denoting the instanton density. From the asymptotic 
behavior and we see immediately, that for the standard zero mode 
Mi{0) > 0, whereas constrained zero modes give Mci{0) = 0. In Fig. 1 we show 
the quark mass function for the constrained and unconstrained zero modes. 
As advocated above, the constrained zero modes lead to (|l|). The momentum 
at which M(p2) starts dropping as p"^ is being decreased is controlled by the 
large-scale parameter R. When R is increased, the constrained profile tends 
to the standard zero mode solution everywhere except the point p^ = 0. 

Let us pass to another approach hinting to Eq. (0), namely QCD sum rules 
with nonlocal condensates Q^x"^) = {q{x)Pexp[ig dzA{z)]q{0)), etc 
The local quark condensates, < qq >, < q{D'^)"-q >, n = 1,2,... appear as 
expansion coefficients of the correlator Q{x'^) in the variable a;2/4. Now, we 
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p2 [GeV2] 

Fig. 1. Dynamical quark mass as a function of the square of the Euclidean momen- 
tum. The dashed line corresponds to the standard zero mode solution of Eq. (^), 
while the solid line corresponds to the constrained zero mode profile, Eq. (^). The 
parameters are p = 1.7GeV~^ and R = 3p. 

make a crucial assumption needed for our argumentation: the quark conden- 
sate and its moments can be identified with the expressions obtained at the 
one-quark-loop level in the effective quark model. This is valid in the \axge-Nc 
limit. We have therefore (for the single-fiavor)^ 



{qq) = -AN,J 



(27r)4p2 + M2(p2)^ 



and 



\ = 



can be found[^ 

.2^ o Mip') 



d'p M{p^ 



{qq)J (27r)4^ p^ + M^{p^)' 



(9) 



(10) 



p'^ + M'^ip'^) 



'111 



where Q{p'^) is proportional to the Fourier transformation of the scalar quark 



2 For simplicity, we are neglecting the effects of the quark wave-function renormal- 
ization by setting Z{q'^) = 1 in Eq. (^). This amounts to neglecting possible vector 
interactions, as is the case of the instanton calculations. 

^ Here and below we silently assume that nonlocal condensates are defined 
in gauge-invariant way by putting in appropriate definitions the path-ordered 
Schwinger phase factor. In QCD sum rules this factor becomes unity throughout 
the use of the Fock-Schwinger gauge. A similar situation occurs in the instanton 
model, where we have to transform the fields from standard (singular or regular) 
gauges to the axial one [13|. This results in modification of Eq. We omit these 
complications in the present work, since the influence of the phase factor is not 
important for the qualitative analysis of the dynamical quark mass. 
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condensate Q{x'^). In QCD sum rules one typically makes a hypothesis for 
Qip^). For instance, one of the ansatze successfully applied in QCD sum rules 

Q(p) = -^^exp(--), (12) 

with average quark virtuality in QCD vacuum being 

a; = 2Al (13) 

Then, the moments of the quark condensate become parameterized in a simple 
way. Thus, in QCD sum rules its the nonlocal quark condensate Q{p^) which 
is the basic quantity, which is in a sense opposite to quark models or the 
instanton approach, where we start with M{p'^). The ansatz (|12D is compatible 
with recent lattice measurements of scalar quark correlator, Q{x^), on the 
lattice PJ^. The dynamical quark mass Mij)^) is related to the nonlocal quark 
condensate Qip^) via Eq. (|TI1). 



We shall now explore the consequences of Eq. (|TTD and try to invert it in order 
to obtain the quark mass as a function of the condensate (assumed to be known 
from the QCD sum rules or lattice simulations 0). A few implications of 
([TT| ) are immediate. We can see that due to the rapid decrease of the profiles 
Qip"^) and M{p'^) at large the relation becomes linear asymptotically: 

M{p^) = ^p'^Qip'^) at large p^. (14) 

This is the dilute instanton medium regime, where the single-instanton effects 
dominate. Next, from (^) it follows that for real valued M{p^) 

\/?Q{p^) < 1 (15) 

for any value of p"^ > 0. Thus, in general, one can consider three different sit- 
uations. First, if the model profile function for the nonlocal condensate Q(p^) 



violates the bound (|T^), then one cannot find a quark- model representation 
for the results in terms of real valued functions. The other two possibilities, 
discussed by us below, correspond to the case when the bound (|15D is satu- 
rated at some points i/p^ = ""^i; or, finally, to the case of the strict inequality 
for all momenta, \fp^Q[j>^^ < 1- Let us first consider the last possibility, where 
we find 



The sign of the root is chosen in such way that the large-p^ asymptotics ([T^) 
is satisfied. The feature of this solution is that M(p2 = 0) = 0. Indeed, at 
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small we get 

M{p') = \q{Q)p' + 0{p'). (17) 

In the case where the bound Eq. (|15|) is saturated at real valued points i/p^ = 
mj, the solution can smoothly inmp from one branch to another. We illustrate 
this for the case of a single point, rriQ. At this point we have the identity 

M{ml) = mo. (18) 

In this case the inverse of (|TTD is 



l-Sgn(p2-m2) Jl-p2Q2(p2) 

M{P ) = . (19) 

The sign function flips the solution from one branch to another at the point 
= rriQ. This prescription ensures the continuity of M{p'^) with all its deriva- 
tives. Note that now, due to branch-switching, we have M(0) = ^l^y > 0. 



We illustrate the above considerations in the model (^21) . The condition ([18|) 
occurs for 



which numerically gives Acrit = 505MeV if one fixes {qq) = — (225MeV)^. At 
A > Acrit we have only one branch. At A = Acrit we have the branch-switching, 
according to ([T9|). The situation is depicted in Fig. 2. 

One determines the values of model parameters, such as A, by fitting observ- 
ables. In quark models one typically fits the pion decay constant, F.^ to its 
experimental value, which in the chiral limit equals 86 MeV |T^. The expres- 
sion for Ft, is P,P 



oo 



Nc f M{uY - uM{u)M\u) + u^M'iuf 



" 47r2 7 {u + M{uf) 

or it may be rewritten in terms of Q as: 



oo 



N, 7^ AQ{uf + AuQ{u)Q'{u)+Au^Q\uf-?,uQ{uf 



^ 647r2 7 1 - AuQiuf 

Above, M'{u) = -^M{u) and Q'{u) = -^Q{u). Figure 3 shows F^ evaluated 
with (|12|), plotted as a function of A. We note that the experimental value of F^, 
favors A ^ 0.55 GeV, comfortably above the critical value, such that indeed 
we have a one-branch situation and, consequently, (p. The average quark 
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p2/A2 

Fig. 2. Dynamical quark mass for the nonlocal condensate given by Eq. ([T^), plotted 
as a function of the square of the Euclidean momentum expressed in units of the 
cut-off parameter A. The curve labeled A > Acrit satisfying condition (|^) corresponds 
to A = 554MeV, which fits the value of the pion decay constant in the chiral limit, 
= 86 MeV. The curve labeled A = Acrit displays the branch-switching. At high 
momenta the solution is given by the solid line, and at low momenta by the dashed 
line. 
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Fig. 3. -Ftt as a function of A (solid line) and the experimental value of 86MeV 
(dashed line). 

virtuality (|T3]) is estimated as = 0.6 GeV^ and fits the value obtained in 
the QCD sum rules = 0.5 ± 0.1 GeV^ [l^. Another interesting feature of 



ansatz ([T^) is that the quark propagator does exhibit analytical confinement 
of quarks. This is because the equation p"^ = M^(— p^) has no real solution at 
positive (Minkowskian) p"^. 



The important feature of effective models is the compliance to symmetries, 
in particular to the chiral symmetry. The consistency of quark models with 
nonlocal regulators with symmetries has been shown in Ref. [|lT| and we are not 
going to repeat it here. Since the proofs never use the explicit form of M(p^), 
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the results are true for any mass function, and thus also hold for functions 
with the property (|l|). 



Finally, we examine the analytic properties of the pion propagator. The pur- 
pose here is to check whether the vanishing quark mass does not induce an 
infrared cut, which would be unphysical. In nonlocal quark models the inverse 
pion propagators is given by |1TT1 , [T^ , [15| , [T6| 



K-\q) = const (23) 



M{p + i)M{p - f ) (p^ - 


3l + M{p + l)M{p-l)) 


M{pf 


((p + ff + M(p+f)2) 


((p - f f + M(p - f )2) 


p2 + M(p)2 



It is easy to verify that the mass factors M(p + |)M(p — |) in the numerator 
cancel the infrared divergence from the denominator, and as a result the pion 
propagators has no cuts at low q. Its analytic structure in the low-energy 
region consists of the pole at = 0, as it should be. If it were not for the 
presence of M{p + ^)M{p — |) in the numerator, expression ( ^3] ) would have 
developed a cut in the variable reaching from to — oo. Analogous result 
holds for the cr-field propagator. 

We conclude that the scenario of a vanishing dynamical quark mass at zero 
quark virtuality is possible, as well as free of pathologies. In forthcoming work 
we are going to study the processes that are sensitive to infrared behaviour of 



nolocal nonperturbative propagator. In particular, in '^Tj it was shown that 
the pion light-cone wave function in the middle region is related to derivative 
of the vacuum nonlocality in the infrared region. Finally, the models discussed 
in this paper lead to, in addition to Eq. analyticity of M(p2) in the vicinity 
of zero. This is not the case of, e.g., standard zero modes, where the function 
develops a cut at p^ = 0, extending along the whole Minkowski region. This 



nonanalyticity becomes an obstacle in soliton calculations in such models |T8 
where other regulators must be used, such that they can be continued to the 
Minkowski momenta near the origin. Models discussed in this paper are free 
of that problem. 

The authors are thankful to A. P. Bakulev, G.V. Efimov, N.I. Kochelev, S.V. 
Mikhailov, S.N. Nedelko, M.K. Volkov for discussions. 
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